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Abstract 

Anomalous quartic couplings between the electroweak gauge bosons may contribute to 
the vv^^ and qq77 final states produced in e''~e~ collisions. This analysis uses the LEP2 
OPAL data sample at centre-of-mass energies up to 209 GeV. Event selections identify 
1/1777 and qq77 events in which the two photons are reconstructed within the detector 
acceptance. The cross-section for the process e'''e~ qq77 is measured. Averaging over 
all energies, the ratio of the observed e"'"e~ qq77 cross-section to the Standard Model 
expectation is 

R(data/SM) = 0.92 ± 0.07 ± 0.04, 

where the errors represent the statistical and systematic uncertainties respectively. The 
Z/I777 and qq77 data are used to constrain possible anomalous W"'"W~77 and ZZ77 cou- 
plings. Combining with previous OPAL results from the W''~W~7 final state, the 95% 
confidence level limits on the anomalous coupling parameters Og, a^, Oq and are 
found to be: 

-0.007 GeV"2 < allt^ < 0.023 GeV"^ 
-0.029 GeV-2 < a^/A^ < 0.029 GeV-^ 
-0.020 GeV-2 < a^/A^ < 0.020 GeV-^ 
-0.052 GeV~2 < a^/A^ < 0.037 GeV-^ 

where A is the energy scale of the new physics. Limits found when allowing two or more 
parameters to vary are also presented. 
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1 Introduction 

In the Standard Model (SM) self-interactions of the vector boson fields arise due to the 
— |W^i,-W'''' term in the electroweak Lagrangian. In addition to the tri- linear couplings, 
this term leads to quartic gauge couplings (QGCs) of the form WWWW, WWZZ, WW77 
and WWZ7. The strength of the coupling at these vertices is specified by the SU(2) x U(l) 
gauge invariant form of the electroweak sector. Studying processes to which these QGCs can 
contribute may therefore yield further confirmation of the non-Abelian structure of the SM or 
signal the presence of new physics at as yet unprobed energy scales. At LEP energies it is only 
possible to probe quartic gauge couplings which produce at most two massive vector bosons in 
the final state. The processes at LEP which are sensitive to possible anomalous quartic gauge 
couplings (AQGCs) are shown in Figure 1. 




Figure 1: The diagrams sensitive to possible anomalous quartic couplings in the e"^e — > 
W"'"W~7, e+e" —>■ uujj and e+e" C[q77 final states. 

The formalism for the extra genuine quartic terms relevant at LEP has been discussed 
widely in the literature [1-7]. Genuine quartic terms refer to those that are not associated with 
any tri-linear couplings, which are already constrained by analyses using the e"'"e~ — > W+W~ 
process. In the parametrisation first introduced in [1] the two lowest dimension terms that give 
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rise to quartic couplings involving at least two photons are: 



CP, = -^Y^F.^^'^W'^''^'^' 

where F'^'^ is the photon field strength tensor. These are C and P conserving and are obtained 
by imposing local U(l)gj^ gauge symmetry, whilst also requiring the global custodial 811(2)^, 
symmetry that preserves the constraint that the electroweak parameter p = 1. We note that 
the custodial SU(2)^ field vector is 



and identifying 




1 

2 cos^ 6w 

yields, in terms of the physical fields, W~ and Zq, 



2 W 2 Z 

p2 „W 2 Z 

Thus, both >Cg and C% generate W"'"W~77 and ZZ77 couphngs, with the parameters and 
now being distinguished for the W and Z vertices to comply with the more general treatment 
in [5]. In all cases the strengths of the quartic couplings are proportional to l/A? where A is 
interpreted as the energy scale of the new physics. 

Limits on AQGCs from LEP data have been published by the OPAL and L3 collaborations 
[8-11]. This paper describes limits on AQGCs obtained by OPAL from the processes e"'"e~ — > 
^'^'77 and e"'"e~ qq77 from all data recorded above the Z pole. For both processes the 
dominant SM background arises from initial-state radiation (ISR). The limits obtained from 
e"'"e~ — > 1/I777 and e"'"e~ qq77 are combined with the limits obtained by OPAL from the 
process e+e~ — > W''"W~7 [11]. 

Since cross-sections for the qq77 final state have not previously been measured explicitly 
by the OPAL collaboration at LEP2, these measurements are presented in this paper and are 
compared with the SM expectation. 



2 The OPAL Detector and Data Samples 

The OPAL detector included a 3.7 m diameter tracking volume within a 0.435 T axial mag- 
netic field. The tracking detectors included a silicon micro- vertex detector, a high precision gas 
vertex detector and a large volume gas jet chamber. The tracking acceptance corresponds to 
approximately | cos^| < 0.95 (for the track quality cuts used in this study) ^. Lying outside the 

-•^The OPAL right-handed coordinate system is defined such that the origin is at the centre of the detector 
and the z axis points along the direction of the e~ beam; 6 is the polar angle with respect to the z axis. 
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solenoid, the electromagnetic calorimeter (ECAL) consisted of 11704 lead glass blocks having 
full acceptance in the range | cos^| < 0.98 and a relative energy resolution of approximately 
6 % for 10 GeV photons. The hadron calorimeter consisted of the magnet return yoke instru- 
mented with streamer tubes. Muon chambers outside the hadronic calorimeter provided muon 
identification in the range | cos6'| < 0.98. A detailed description of the OPAL detector can be 
found in [12]. 

From 1995 to 2000 the LEP centre-of-mass energy was increased in several steps from 130 to 
209 GeV. For the analysis of the qq77 channel, this entire data sample is used, corresponding 
to 712 pb^\ The 1/7777 analysis is restricted to 652 pb^^ of data recorded above 180 GeV. 
The integrated luminosities at each centre-of-mass energy for the vV^^ analysis are lower than 
those for the qq77 analysis due to tighter requirements on the operational status of the detector 
components. 

3 Monte Carlo Models 

A number of Monte Carlo (MC) samples, all including a full simulation [13] of the OPAL 
detector, are used to simulate the SM signal and background processes. For the vV^^ final 
state NUNUGPV [14] is used to model both the dominant SM doubly-radiative return process 
and the supplementary AQGC processes, with KK2F [15] being used as a cross-check on the SM 
expectations. For the qq77 final state, the KK2F program is also used. For the background 
processes, the concurrent MC tandem [16] of KORALW and YFSWW is used to simulate 
the background from four-fermion final states with fermion fiavour consistent with being from 
W+W" final states. The KORALW program [17] is used to simulate the background from 
four-fermion final states which are incompatible with coming from the decays of two W-bosons 
{e.g. e+e~ qq^u+yU"). For both signal and background processes JETSET [18] is used to 
model the fragmentation and hadronisation of final state quarks. The two-fermion background 
process e"'"e~ — > Z/7 — > t^t~ is simulated using KK2F. The background in the qq77 event 
selection from multi-peripheral two-photon diagrams is negligible. The WRAP program [7] is 
used to determine the effects of AQGCs in the qq77 channel. 

4 The 1/I777 Final State 

4.1 vV^^ Event Selection 

The selection proceeds in two stages: 

Acoplanar photon pair selection: This event selection employs standard criteria described 
in detail elsewhere [19,20]. Candidate events must meet the kinematic requirement of there 
being at least two photons, either both with energy E.^ > 0.05£'beam and polar angle 9j satisfying 
I cos 6*^1 < 0.966, or one with Ey > 0.05£Jbeam, | cos^^| < 0.966 accompanied by a second with 
E.y > 1.75 GeV, |cos^^^| < 0.8 that has an associated in-time time-of-flight detector signal. 
Events with three final state photons (e"'"e~ —>■ 1/17777) are permitted, the subsequent selection 
criteria then being applied to the two photons with the highest reconstructed energies. The 
system consisting of the two highest energy photons must have a momentum transverse to the 
beam axis, , satisfying > 0.05£^beam- Additional requirements are then made on the 
photon conversion consistency (charged track veto), the electromagnetic calorimeter cluster 
shape, the forward energy vetoes and the muon vetoes. The e"*'e~ —>■ 77(7) background is 
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suppressed whilst retaining the events with missing energy by imposing further cuts on the 
energies and angles of the selected two or three photon system. These include the requirements 
that the total energy in the electromagnetic calorimeter does not exceed 0.95y^ and also that 
the acoplanarity^ angle of the two highest energy photons be greater than 2.5°. 

The efficiency for SM e"'"e~ — > 1/1/77(7) events within the kinematic acceptance of the acopla- 
nar photon pair selection is approximately 66% [20]. The expected background contribution 
from processes other than e+e" — > 1/1/77(7) is less than 1% [19,20]. 

Suppression of Standeird Model background: To suppress the SM contribution, princi- 
pally the forward-peaked doubly-radiative return process, the following additional cuts are 
applied to the events passing the acoplanar photon pair selection: 

• The two highest reconstructed photon energies, E^i and E^2, niust both be greater than 
10 GeV. This cut has little effect on any AQGC contribution, which gives rise predom- 
inantly to photons of high energy, but does suppress the doubly-radiative return back- 
ground. 

• |cos^-yi| < 0.9, |cos^^2| < 0.9, where again the subscripts refer to the two photons with 
highest reconstructed energy. This requirement further suppresses the doubly-radiative 
return background, which is forward peaked as expected for initial-state radiation pho- 
tons. 

These cuts were optimised on SM MC to yield the maximum sensitivity to the anomalous 
couplings. 

4.2 Sensitivity of e+e" — ^ vV")") to Anomalous QGCs 

Table 1 lists the number of data events accepted by the Z/Z777 event selection compared to the 
SM expectation, binned by centre-of-mass energy. There is excellent agreement between the 
predictions of NUNUGPV and the KK2F MC program [15] used as a cross-check. The SM 
predictions describe the data well. 

Approximately 4.0—4.7% of real data events, depending on the centre-of-mass energy, are 
expected to fail the acoplanar selection due to the effects of random coincidental activity. These 
rates have been evaluated from samples of random beam-crossing events collected throughout 
the data-taking periods. All quoted MC accepted cross-sections have been corrected for these 
unmodelled effects. 

For the selected events. Figure 2 shows the distribution of the invariant mass recoiling 
against the photons, Mrcc, and the distribution of the energy of the photon with the second 
highest reconstructed energy, E^2- In both cases the data arc well described by the SM ex- 
pectation. Figure 2 also shows the effects of anomalous couplings on these distributions. For 
the recoil mass, increasing the coupling at the ZZ77 vertex increases the cross-section at the Z 
mass peak, whereas the effect of the W^W~77 vertex can mainly be seen in the low recoil mass 
region of the plot. Similarly, the effects of the different quartic vertices can be distinguished in 
different regions of the E^2 distribution. 

Constraints on AQGCs are derived employing a maximum likelihood fit that uses bins in the 
Mj-ec and E^2 distributions at each centre-of-mass energy. The ten bins are defined in Table 2, 

^The acoplanarity angle is defined as tt minus the opening angle between the two photons when projected 
onto a plane perpendicular to the beam axis. 
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jCdt 


Data 


SM Expectation 


[GeV] 


[pb-i] 




NUNUGPV 


KK2F 


180-185 


53.9 





2.5 


2.5 


188-190 


175.2 


10 


7.9 


7.9 


191-192 


28.8 


1 


1.3 


1.3 


195-196 


71.6 





3.1 


3.0 


199-201 


73.7 


3 


3.0 


2.9 


201-203 


36.7 


1 


1.5 


1.4 


203-209 


210.6 


5 


8.3 


8.0 


Total 


652 


20 


27.6 


27.0 



Table 1: Numbers of i/V^^ events passing the event selection by centre-of-mass energy compared 
to the SM expectations from both KK2F and NUNUGPV. All MC accepted cross-sections have 
been corrected for efficiency losses due to random coincident detector hits. 



together with the corresponding numbers of events observed and expected in the SM summed 
over centre-of-mass energies. The choice of binning reflects the differing effects of the anomalous 
couphngs on the different regions of the M^ec and E^2 distributions and was optimised on SM 
MC for maximum sensitivity to the coupling parameters, inclusive of systematic effects. 



Bin Number 




[GeV] 


E^2 [GeV] 


Observed 


Expected 


1 


< 


60 


10-25 





0.1 


2 


< 


60 


25-45 





<0.1 


3 


< 


60 


>45 





<0.1 


4 


60- 


-80 


10-25 


1 


0.5 


5 


60- 


-80 


25-45 


2 


0.4 


6 


60- 


-80 


>45 





0.1 


7 


80- 


- 120 


10-25 


5 


11.7 


8 


80- 


- 120 


25-45 


6 


8.3 


9 


80- 


- 120 


>45 


1 


0.8 


10 


> 


120 


>10 


5 


5.7 


Total 


20 


27.6 



Table 2: The binning of the likelihood function for the z/z/77 events together with the corre- 
sponding numbers of events observed and expected in the SM. 

4.2.1 Systematic Uncertainties (1^1777) 

The systematic errors in this analysis are found to be small in comparison to the statistical 
error from the 20 selected data events. 

Experimental uncertainties: The main experimental systematic uncertainty arises from the 
accuracy of the modelling of the energy scale and resolution of the electromagnetic calorimeter. 
The evaluation of this is based on a comparison of reconstructed events with two beam-energy 
photons in the final state e+e~ 77 with those simulated in MC. Additional degradations 
in the resolution and scaling were then applied to the accepted SM cross-sections (both total 
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and in the analysis bins) to evaluate the systematic uncertainties, separately for the barrel 
(I cos^^^l < 0.7) and end-cap (0.7 < | cos^^^| < 0.9) regions of the detector and for each year 
of data taking. These uncertainties result in relatively large fractional systematic uncertainties 
for individual analysis bins (approximately 20% for the bins with smallest cross-section, i.e. 
bins 2 and 3 of Table 2) though these propagate through to small overall errors of less than 
1 % on the total cross-sections. Possible biases in the measured photon angle were found to be 
neghgible. 

Theory shape uncertainty: The shapes of the SM M^-^c and E^2 distributions from KK2F 
and NUNUGPV have been compared in order to evaluate any possible theoretical uncertainty 
in the SM prediction. Again, the variations in the total cross-sections were small (< 4%), 
but large fractional variations could be seen for bins 1—3 which were hardly populated by the 
statistics available from KK2F. 

Normalisation uncertainty: Other sources of systematic uncertainty have been considered 
and affect primarily the overall normalisation. The uncertainty related to the modelling of 
initial-state radiation (ISR) has been assessed by turning off ISR with finite p^j leading to 
a ±5 % normalisation uncertainty. The cross-sections for NUNUGPV have been compared 
with the predictions of Belanger et al. [5] and the difference used to estimate a normalisation 
systematic uncertainty of ±4 %. In addition, the luminosity error is ±0.3 %. These errors are 
added in quadrature to give an estimate of the overall normalisation uncertainty of 6.4 % which 
is taken to be independent of energy. 

At all centre-of-mass energies and for any combination of the couplings, the available 
NUNUGPV MC statistics amounts to at least one thousand times the data statistics and 
the related MC statistical error is negligible. Similarly, due to the large sample sizes of random 
events analysed, the uncertainties on the corrections for losses due to coincidental random detec- 
tor hits are less than 1 % and are neglected. The systematic error associated with the expected 
background contribution from processes other than e''"e~ — > 1/1/77(7) is also negligible. 

4.3 Limits on Anomalous QGCs from e"'"e~ — > vV^^ 

At each centre-of-mass energy, 15 samples of 2 000 events with differing values of 
and have been simulated. The extra Lagrangian terms are linear in the anomalous cou- 
plings. Consequently, the cross-section has a quadratic dependence and these 15 samples are 
sufficient to parametrise fully (T(a^, Oq, a^, a^). The generated events arc rcweighted using 
matrix element weights from NUNUGPV to obtain Monte Carlo samples corresponding to any 
combination of the anomalous QGCs (a^, Oq, a^, a^). 

For the vV^^ final state, fits for each of the AQGC parameters have been performed to the 
data by summing the likelihood curves obtained from the seven centre-of-mass energies con- 
sidered. The effects of systematic uncertainties are included in the fits. The fitted AQGCs are 
compatible with zero and the resulting 95% confidence level (C.L.) intervals on the anomalous 
couphngs varied individually are listed in Table 3. These limits include the effects of systematic 
uncertainties. The corresponding likehhood curves are shown in Figures 3a-3d. The results of 
a fit allowing two AQGC parameters to vary simultaneously are shown in Figure 4, again with 
the two parameters not plotted fixed at zero. Since anomalous ZZ77 and W"'"W~77 couplings 
affect different regions of the invariant mass and second photon energy distributions, the limits 
on and Oq are largely uncorrelated. The same is true for the limits on and . 
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95 % C.L. Limit 


7/7 /O/O/ 

1/1/ y y 




— U.UUy Vort; V 


"2 ^ / A 2 ^ n nofi PoV" 


2 


777 yo/o/ 




— U.Uo4 lorC V 


< d^j-l^ < U-UoyLreV 


2 






— n 040 OpV" 

U .UtiU VJC V 


-2 , W / A 2 ^ n no7 rtpV" 


-2 




w 

(7 


-0.114 GeV- 


■2 < aW/A2 < 0.103 GeV" 


-2 


qq77 


«0 


-0.012GeV" 


"2 < a^/A2 < 0.027 GeV" 


2 


qq77 




-0.036 GeV 


< a^/A^ < 0.034 GeV- 


2 


W+W-7 


"0 


-0.020 Ge 


< aW/A2 < 0.020 GeV" 


-2 


W+W-7 


w 


-0.053 Ge V- 


■2 < aW/A2 < 0.037 GeV" 


-2 



Table 3: The 95 % C.L. limits on the anomalous QGCs from the OPAL LEP2 data from the 
processes shown in Figure 1. The 1/1/77 ^^'^ QQTT results are described in this paper. The 
limits from the process e"'"e~ W"'"W~7 are described in Reference [11]. All limits include 
systematic uncertainties and correspond to the case where only the coupling in question is 
varied from zero. 

5 The qq77 Final State 

In the SM, photons in the process e"'"e~ — > qq77 are radiated from cither the initial or final 
state fermions. Photons from ISR tend to be produced along the beam direction. Photons 
from final state radiation (FSR) tend to be produced almost coUinear with the quarks and are 
often lost within hadronic jets. For the measurement of the qq77 cross-section a theoretical 
acceptance is defined which is well matched to the experimental sensitivity. The cross-section 
is defined within a qq invariant mass region dominated by the Z exchange diagrams. 

The e"''e~ — > qq77 cross-section measured in this paper corresponds to the following accep- 
tance with respect to the qq77 system: 

• There must be at least two photons satisfying: 

I) E^i > 5 GeV, where E^i is the energy of photon i, 

ii) |cos^-yj| < 0.95, where O^i is the polar angle of photon i, 

in) cos < 0.90, where 91^^ is the angle between photon i and the direction of the 
nearest quark. 

• |Mqq - Mz| < 3rz. 

The quantity Mqq is defined as the propagator mass, i. e. the invariant mass of the qq system 
before FSR. Photons from FSR are not considered as signal and interference between ISR and 
FSR is neglected. 

5.1 qq77 Event Selection 

The selection of the qq77 events proceeds in three stages: 

e+e~ — >■ qq event selection: e+e" —>■ qq events are selected using the algorithm described 
in [21]. 
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Photon identification: Photon candidates can be identified as either unassociated electro- 
magnetic calorimeter (ECAL) clusters or photon conversions, following the procedure described 
in [11]. Only photons with measured energy > 5 GeV and polar angle |cos^^| < 0.95 are 
retained. The remainder of the event is forced into two jets using the Durham algorithm [22]. 
Finally, to reduce background from photons from the decays of neutral hadrons, e.g. 7r° and 
T] decays, the photons are required to be isolated from the reconstructed jets by requiring 
cos ^^-jET < 0.9, where ^^-jet is the angle between the photon and the direction of the closest 
reconstructed jet. Photon candidates which fail this isolation criterion are merged to the nearest 
jet and the jet energy is recalculated. Events with two or more identified photons satisfying the 
above requirements are retained for the analysis. For photons within the MC generator level 
acceptance > 5 GeV, |cos^^| < 0.95 and cos^^q < 0.9, the photon identification efficiency 
is about 88%. The requirement of two identified photons therefore rejects approximately 23% 
of the qq77 signal. 

Kinematic requirements: The reconstructed mass of the hadronic system, Mqq, is required 
to be consistent with Mz- For about 90 % of the events Mqq is obtained from a kinematic fit 
which imposes the constraints of energy and momentum conservation. In the first instance the 
fit assumes a four-body final state consisting of two jets and two photons. If the fit probability 
is less than 0.01, the fit is performed allowing for an unobserved photon along the e+e" beam 
axis. For events where this fit probability is also less than 0.01, the hadronic mass is taken 
to be the recoil mass calculated from the reconstructed momenta of the two photons. The 
number of events with mass reconstructed in the three possible categories is consistent with 
MC expectation. The reconstructed invariant mass spectrum before the cut on Mqq is shown in 
Figure 5. Events within the region 75 GeV < Mqq < 125 GeV are considered qq77 candidates. 
The cut on Mqq removes 47 events in the data compared to the SM expectation of 58.6. Due to 
experimental resolution this mass window is larger than that used in the kinematic definition of 
the cross-section. Nevertheless, this cut rejects approximately 6 % of the qq77 events satisfying 
the signal definition. 

After applying the cut on Mqq a total of 176 events are identified in the data, consistent 
with the SM expectation of 191.0. Figures 6a-6e show the distributions of Eji, E^2, |cos^-yi|, 
|cos^^2| and E^i + E^2 for selected events. Figure 6f shows the distribution of the maximum 
I cos I of the two highest energy photons in the event. In each case the data are in good 
agreement with the SM expectation. 



5.2 Cross-section Results 

The qq77 cross-section is determined within the above acceptance definition. Cross-section 
values are obtained for the seven different centre-of-mass energy ranges listed in Table 4. The 
qq77 cross-section is calculated from 

C'qq77 — r - T ' 

where N^hs is the accepted number of events, A|^^ is the SM expected number of background 
events, and L — J Cdt is the integrated luminosity, given in Table 4. The qq77 selection 
efficiency, Sqq^^, is evaluated using the KK2F MC samples and includes feed-through from 
genuine qq77 events outside the signal acceptance (a contribution of approximately 12%). 

The numbers of events selected at each energy are listed in Table 4 along with the quantities 
used to calculate the cross-sections. Also shown are the derived cross-sections for the above 
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signal acceptance. The systematic uncertainties are described below. The results are consistent 
with the SM expectation, as shown in Figure 7. Averaging over all energies, and taking into 
account correlated systematic uncertainties the ratio of the observed to expected cross-sections 
is 

i?(data/SM) = 0.92 ± 0.07 ± 0.04, 



where the errors represent the statistical and systematic uncertainties respectively. 





< > 








iVobs 




aqq^^(SM) 


[GeV] 


[GeV] 


[pb-i] 


[%] 






[fb] 


[fb] 


130.0-137.0 


133.0 


10.6 


76.2 ±4.0 


1.1 ±0.3 


8 


848 ± 350 ± 57 


738 


160.0-173.0 


166.9 


20.3 


79.4 ± 3.2 


1.0 ±0.2 


5 


247 ± 139 ± 17 


412 


180.0-185.0 


182.7 


57.2 


77.5 ± 3.1 


2.7 ±0.5 


10 


164 ± 71 ± 13 


333 


188.0-189.0 


188.6 


183.1 


77.7 ±2.9 


9.5 ±1.6 


53 


305 ±51 ± 16 


309 


191.0-196.0 


194.4 


105.7 


77.4 ±2.9 


4.3 ±0.7 


25 


254 ±61 ± 13 


288 


199.0-204.0 


200.2 


114.1 


78.4 ±2.9 


3.0 ±0.6 


26 


257 ± 57 ± 12 


270 


204.0-209.0 


205.9 


220.6 


76.0 ±2.9 


7.2 ±1.3 


49 


250 ± 47 ± 12 


257 



Table 4: Selected qq77 events and cross-section results for the seven different y/s ranges used in 
the analysis. The y/s range, the mean luminosity weighted value of y/s and the corresponding 
integrated luminosity, / Cdt, are listed. For the measured cross-sections, the uncertainties are 

respectively statistical and systematic. The uncertainties on the efficiencies and backgrounds 
are the estimated systematic uncertainties including a contribution from finite MC statistics. 
Also shown is the SM expectation from KK2F. 

5.2.1 Systematic Uncertainties (qq77) 

The systematic uncertainties on the qq77 selection efficiency and on the expected number 
of background events are estimated to be 2.7% and approximately 20% respectively. The 
systematic uncertainties, described below, were obtained in the same manner as described in 
Ref. [11] where further details may be found. In addition the contributions to the systematic 
uncertainties due to finite MC statistics are included in the numbers listed in Table 4. 

Photon identification and isolation: A systematic uncertainty of 1 % is assigned to cover 
the uncertainties in the simulation of the photon conversion rate and the accuracy of the 
simulation of the electromagnetic cluster shape [23]. The systematic error associated with the 
isolation requirements depends on the accuracy of the MC simulation of the fragmentation 
process in hadronic jets. This is verified in Z — > qq events recorded at y/s ~ Mz during 
1998 — 2000. For each selected event, the inefficiency of the isolation requirements is determined 
for cones of varying half-angle defined around randomly orientated directions. The inefficiency 
of the isolation cuts is parametrised as a function of the angle between the cone and the nearest 
jet. For all cone half-angles the inefficiency in the MC and data agree to better than 1 %; 
consequently a 1 % systematic error is assigned. These two effects give a total uncertainty on 
the identification efficiency for a single photon of 1.4 %. Since two photons arc required in the 
analysis of qq77 this corresponds to an uncertainty in the qq77 efficiency of 2.8%. 
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Photon energy scale and resolution: A bias in the energy scale for photons (data relative 

to MC) in the region of the energy cut, i.e. Ey 5 GeV, would result in a systematic bias in 
the qq77 cross-section measurement. The uncertainty on the ECAL energy scale for photons in 
this region is estimated by examining the invariant mass distribution of pairs of photons from 
7r° decays in e+e~ — > qq events recorded at ^/s ~ Mz during 1998—2000 and e+e~ qq(7) 
events recorded at a/s > 180 GeV. As a result a 4% systematic uncertainty on the ECAL 
energy scale in the region of -E-y ~ 5 GeV is assigned. The resulting systematic uncertainty on 
the qq77 cross-section is 1.5%. 

The systematic error from the uncertainty in the ECAL energy resolution is obtained in 
a similar manner to that used for the ECAL energy scale using the same 7r° sample. There 
is no evidence for a statistically significant difference between the energy scales in data and 
MC. The statistical precision of the comparison, ±10 %, is used to assign the energy resolution 
uncertainty which, when propagated to the uncertainty on the qq77 cross-section, yields a 
systematic error of ±0.6 %. 

Photon angular acceptance: The systematic error associated with the acceptance require- 
ment of |cos 6y \ < 0.95 depends on the accuracy of the MC simulation of the angular reconstruc- 
tion of ECAL clusters at the edge of the acceptance. By comparing the reconstructed polar 
angle of leptons from different detectors (ECAL, tracking, muon chambers) in e+e~ — > e+e~ 
and e'^e" — > //"^/x" events the ECAL acceptance is known to ±3 mrad. This uncertainty results 
in a 0.6% uncertainty in the qq77 cross-section. 

Background uncertainties (AT^^): The dominant source of background is from e"'"e~ — > 
Z/7 ^ qq7 where one of the identified photons is from ISR and the other is associated with the 
hadronic jets. A photon associated with the hadronic jets may be either from FSR in the parton 

shower or from the decay of a hadron {e.g. n or rj decays). From the studies presented in [11] 
a 30 % systematic uncertainty on this background contribution is assumed. The systematic 
uncertainties on the small background contributions from four-fermion events and from tau- 
pair events are negligible. An additional 0.8% error is assigned to cover uncertainties in the 
e+e~ — > qq selection. 

5.3 Limits on Anomalous QGCs from e+e~ qq77 

The e"^e~ — > qq77 process is sensitive to the anomalous ZZ77 vertex and the possible couplings 
Cq , a^. To set limits on these a binned maximum likelihood fit to the observed distribution of Ej2 
is performed in 5 GeV bins. Fits are performed to the data for the seven separate energy ranges 
of Table 4 and the resulting likelihood curves are summed. The effects of anomalous couplings 
are introduced by reweighting events generated with KK2F using the ratio of anomalous QGC 
to SM matrix elements obtained from the WRAP program [7] . The resulting likelihood curves 
for one-dimensional fits to Og and separately are shown in Figures 3a and 3b. From these 
curves, 95 % C.L. upper limits on the anomalous couplings are obtained, shown in Table 3. 
The limits include the effect of the experimental systematic errors and assume a 5 % theoretical 
uncertainty (obtained by comparing the predictions of KK2F and WRAP over the centre-of- 
mass range considered in this pubhcation). The 95 % C.L. contour obtained from a simultaneous 
fit to Oq and is shown in Figure 4a. 
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6 Combined Limits on Anomalous QGCs from the qq77, 
vT>^^ and W"'"W~7 processes 



The summed one-dimcnsional likelihood curves for the parameters Oq and from the qq77 
and vV^^ final states are shown in Figures 3a and 3b. In this combination the small effect 
of correlated systematic uncertainties between the two channels has been neglected^. The 
corresponding combined 95 % confidence level limits on possible anomalous contributions to 
the ZZ77 vertex are 

-0.007 GeV-2 < allk^ < 0.023 GeV"^ 
-0.029 GeV-2 < a^/A^ < 0.029 GeV'^ 

When both ZZ77 parameters are allowed to vary simultaneously the likelihood contours of 
Figure 4a are obtained. 

The limits on possible anomalous contributions to the WW77 vertex obtained here from 
the Z/I777 channel are combined with the previous OPAL limits from the e"^e~ — > W+W^7 
process [11]. The resulting likelihood curves are shown in Figures 3c and 3d, again assuming 
the systematic uncertainties for the two channels are uncorrelatcd. The corresponding 95% 
confidence level limits on anomalous contributions to the W"'"W~77 vertex are: 

-0.020 GeV-2 < aW/A2 < 0.020 GeV-^ 
-0.052 GeV-^ < a^/A.'^ < 0.037 GeV"^. 

The likelihood contours for these two parameters are shown in Figure 4b. 

In the hterature the assumption that of = has been made (see for example Ref . [4] ) . The 
validity of the hnking of the W''"W~77 and ZZ77 couplings has been questioned in Ref. [5]. For 
completeness, limits are presented for the case where af = of by combining the one-dimcnsional 
likelihood curves from the z/1777, qq77 and W"'"W~7 processes, shown in Figures 3e and 3f. 
The combined likelihood yields the 95 % confidence level limits: 

+0.002 GeV-2< a^/A^ < 0.019 GeV-^ 
-0.022 GeV-^ < a^/A^ < 0.029 GeV"^. 

The corresponding two-dimensional fit is shown in Figure 8. 

7 Conclusion 

Event selections for the processes 1/I777 and qq77 are presented. The selected qq77 events are 
used to measure the cross-section for the process e"'"e~ — > qq77. Averaging over all energies, 
the ratio of the observed e+e" — > qq77 cross-section to the Standard Model expectation is 

R(data/SM) = 0.92 ± 0.07 ± 0.04, 

where the errors represent the statistical and systematic uncertainties respectively. The selected 
1/1/77 and qq77 events are used to constrain possible anomalous W^W^77 and ZZ77 couplings. 

^The correlated component of the systematic uncertainty on the event selection efficiencies is estimated to 
be 2 %, dominated by correlated uncertainties from the photon energy scale and photon angular acceptance. 



13 



When these results are combined with previous OPAL results from the W+W^7 final state the 
95 % confidence level limits on the anomalous coupling parameters Oq , a^, and are found 
to be: 

-0.007 GeV-^ < ag/A^ < 0.023 GeV-^ 
-0.029 GeV-2 < a^/A^ < 0.029 GeV-^ 
-0.020 GeV-^ < a^/A2 < 0.020 GeV-^ 
-0.052 GeV-2 < af/A^ < 0.037 GeV-^ 

where A is the energy scale of the new physics. Limits allowing two or more parameters to vary 
are also presented. 
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Figure 2: Distributions of a) M-^^c and b) E^2 for the accepted Z/I777 events. The points show 
the 180 — 209 GeV data and the histograms show the MC expectation. The hatched histogram 
represents the SM scenario whilst the expected distributions for possible ZZ77 and W'*"W~77 
AQGC hypotheses are shown by the dashed and dotted lines, respectively. 
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Figure 3: Plots a) and b) show the one dimensional minus log hkelihood curves for Oq and 
from the pV'yy channel (dotted line), the qq77 channel (dashed line), and the two channels 
combined (continuous line). Plots c) and d) show the one dimensional likelihood curves for 
and from the uujj channel (dotted line), the W+W~7 channel (dashed line), and the 
two channels combined (continuous line). Figures e) and f) show the combined limits assuming 
Qq — and — . e) The one dimensional likelihood curve for — aQ — (continuous 
line) with the contribution from the Qq from the qq77 and uV'j'y channels (dotted line) and from 
the limit on from the W^W~7 and uV'-f-f channels (dashed line), f) The one dimensional 
likelihood curve for = = (continuous line) with the contribution from from the 
qq77 and uujj channels (dotted line) and from the limit on from the W"'"W~7 and i/V^^ 
channels (dashed hue) . All hkelihood curves include the effects of systematic uncertainties and 
correspond to the case where only the coupling in question is varied from zero. 
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Figure 4: a) The 95% confidence region in (ag,a^) from tlie uV'j'y cliannel (dotted line), tlie 
qq77 cliannel (dashed line), and the two channels combined (continuous line), b) The 95% 
confidence region in {a^,a^) from the 1/1/77 channel (dotted line), the W"'"W~7 channel [11] 
(dashed line), and the two channels combined (continuous line). In b) the limits from the 
W"'"W~7 channel dominate to such an extent that the limits from the W+W~7 channel alone 
almost coincide with the combined limit. In both a) and b) the position of the best fit (minimum 
of the — InL surface) is indicated by the star and the SM expectation at (0, 0) is shown by the 
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Figure 5: Invariant mass of the hadronic system, Mqq, in selected qq77 events. The arrows 
indicate the cuts used to select the final qq77 sample. The singly hatched histogram indicates 
the background from qq7 events and the doubly hatched histogram (barely visible) indicates 
the small four-fermion and tau-pair backgrounds. 
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Figure 6: Distributions of E^i, E^2, |cos^^^i|, | cos 6*^2!, -£^71 + and the maximum of | cos 6*^1! 
and I cos 6*^2 1 for selected qq77 events. The points show the 130 — 209 GeV data and the 
histograms show the MC expectation. The singly hatched histogram indicates the background 
from qq7 events and the doubly hatched histogram indicates the four-fermion and tau-pair 
backgrounds. The expected distributions for an anomalous QGC parametrised by Oq/A^ = 
0.015 GeV~^ are shown by the dotted lines. 
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